Gravitational instabilities of superspinars 
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that these objects are unstable independently of the boundary conditions, and that the instability 
is strongest for relatively small values of the spin. Also, we give a physical interpretation of the 
various instabilities that we find. Our results (together with the well-known fact that accretion tends 
to spin superspinars down) imply that superspinars are very unlikely astrophysical alternatives to 
black holes. 
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range. Therefore, the observation of rapidly spinning ul- 
tracompact objects could potentially reveal or rule out 
the existence of superspinars. 

One argument against the existence of superspinars 
was put forward in Rcf. 0- There, the authors con- 
structed a toy model for a superspinar by assuming that 
the external surface of the superspinar can be modeled 
as a perfect mirror, i.e. that the reflection coeflicient 
7?. = 1 for waves incident on the superspinar. In this 
case superspinars are destabilized by superradiant ef- 
fects, i.e. by the ergoregion instability first discussed 
by Friedman, Schutz and Comins Q. The ergore- 
gion instability occurs on a dynamical timescale, pos- 
ing a serious challenge to the existence of these objects 
in nature. However, a perfectly reflecting surface may 
be an unrealistic assumption. In general we would ex- 
pect a frequency-dependent reflection coefficient Tl{uj), 
and correspondingly a frequency-dependent transmission 
coefficient T{uj) = 1 — TZ{uj). The exact form of TZ{uj) 
depends on the specific model, but unfortunately no 
exact solutions describing four-dimensional superspinars 
are known. 

A different instability was recently discussed by Dotti 
et al. [m ■ These authors studied perturbations of 
a Kerr solution with a/M > 1 (i.e., unlike Ref. 0, they 
considered an actual naked singularity). They cast the 
linearized perturbation equations in the form of a self- 
adjoint operator and analyzed the discrete spectrum of 



Superspinars are vacuum solutions of the gravitational 
field equations whose mass M and angular momentum 
J = aM violate the Kerr bound, i.e. a > M (here 
and elsewhere in this paper we use geometrical units: 
G = c = 1). These geometries could result from high- 
energy corrections to Einstein's theory of gravity, such 
as those that would be present in string-inspired models 
String-inspired corrections may require a modifica- 
tion of the metric (or some sort of "excision" ) in a small 
region surrounding the curvature singularity at the ori- 
gin, in such a way as to "dress" the singularity. While 
stable stars with a > M are in principle allowed in gen- 
eral relativity^, superspinars have been proposed as an 
alternative to black holes (BHs), and they are therefore 
imagined to have a compactness comparable to that of 
extremal rotating Kerr BHs and to exist in any mass 
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^ Typical equations of state usually lead to stars with a/M < 0.7 
2, that can be treated within a slow-rotation approximation 
3] . However, stable, differentially rotating polytropic stars with 
a/M Ri 1.1 can be produced (e.g.) with the Whisky code @, 
Also, note that the Kerr bound can be easily violated by non- 
compact objects such as the Earth (M/R ~ 7 X lO^^''), which 
has J/M2 ~ 10^. 
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this operator, proving the existence of an infinite number 
of unstable modes . 

Here we generahze the stabihty analyses of Refs. 
and [llj focusing on a superspinar model obtained by 
considering the Kerr solution with a > M. Besides ex- 
tending the study of Ref. 0, we also impose an alter- 
native (and perhaps more physical) prescription for the 
external surface of a four-dimensional superspinar. We 
assume that a perfectly absorbing surface (a "stringy 
horizon") is created by high-energy effects at some ra- 
dius r = tq, and we impose that the reflection coeffi- 
cient Tl{uj) = at that radius. These purely ingoing 
boundary conditions at r = ro are designed to make su- 
perspinars as stable as possible against the ergoregion 
instability of Ref. This instability occurs because, 
when the boundary at r ~ j'q is purely reflecting, the 
negative-energy modes which exist in the ergoregion can 
only leak to spatial infinity by tunneling through a po- 
tential barrier. Modes propagating outside the ergore- 
gion have positive energies. This results in the negative 
energy of the ergoregion modes to decrease indefinitely, 
so that their amplitude becomes unbound, triggering an 
instability. By imposing purely ingoing boundary condi- 
tions at r = ro, we basically allow the negative energy 
trapped in the ergoregion to "fall down a sink" ; this could 
quench the instability to some extent. A similar quench- 
ing occurs for Kerr BHs with a < M, the role of the sink 
being played by the BH horizon. Clearly, if the refiection 
coefficient < 7?,(a;) < 1 the quenching would be less 
efficient. Therefore we conjecture that if supcrspinars 
are unstable when TZ{u}) = 0, it should be impossible to 
stabilize them using any other choice of boundary condi- 
tions. 

In this paper we analyze the stability of supcrspinars 
by imposing either perfectly absorbing (7?.(w) = 0) or 
perfectly reflecting (7?.(w) = 1) boundary conditions at 
some arbitrary radius r ~ tq. We find that, quite inde- 
pendently of ro and of the chosen boundary conditions, 
supcrspinars are unstable to linearized gravitational per- 



turbations. For purely ingoing boundary conditions the 
instability is slightly weaker than in the perfectly refiect- 
ing case, but it still occurs on a dynamical timescale 
T ~ M, i.e. T 5 X 10"^ s for an object with M = Mq 
and T ^ 5 s for a supermassive object with M ^ lO^M©. 
We also show that this result is valid for a wide class 
of theories of gravity. Our findings undermine several 
claims made in the literature that supcrspinars might be 
detected because the shadow they cast due to gravita- 
tional lensing [l^, [I!] or their accretion properties [11, 
[Tsl are different from Kerr BHs with a < M. While this is 
true, supcrspinars are plagued by multiple gravitational 
instabilities, and therefore they are unlikely to be astro- 
physically viable BH candidates. 
I. A SIMPLE MODEL OF SUPERSPINAR IN 
FOUR DIMENSIONS 

Following Gimon and Horava we model a super- 
spinar of mass M and angular momentum J = aM by 
the Kerr geometry 
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where S = r^ + cos^ 6 and A = r^ + - 2il/r. Unlike 
Kerr BHs, supcrspinars have a/A/ > f and no horizon. 
Since the domain of interest is — oo < r < -|-oo, the 
spacetime possesses naked singularities and closed time- 
like curves in regions where g^^ < (see e.g. [Toj). 

We study linear perturbations around the Kerr metric 
(H]) . Using the Kinnersley tetrad and Boyer-Lindquist co- 
ordinates, it is possible to separate the angular and radial 
variables [20| . Small perturbations of a spin-s field arc 
then reduced to the radial and angular master equations 
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where x = cos 6, K = {r'^ + a'^)uj — am and the separation tively. The oscillation frequencies of the modes can be 



constants A and s^Zm are related by 



A = sAim + c^io^ — 2amuj . 



found from the canonical form of Eq. ([2]). Switching 
to a "tortoise coordinate" defined by the condition 
(4) dr^/dr = [r^ + a^)//\, we get 



The equations above describe scalar, electromagnetic and 
gravitational perturbations when s = 0, ±1, ±2 respec- 



(5) 
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where 

^ _ - 2is{r - M)K + A{4irujs - A) dG 

and G = s{r~M)/{r'^ + a'^)+r/\{r'^ + a^)-'^. The eigen- 
values sAim in Eq. ^ can be expanded in a power series 
in the parameter auj as [2l| 

The absence of ingoing waves at infinity imphes [2^ 

r - r-^e'"''- , r-^oo. (7) 

The boundary conditions at r = rg are crucial. Ref. 
assumed a perfect mirror at r = ro, i.e. F(ro) = 0. If in- 
stead we assume the existence of some "stringy horizon" 
at rp, we must impose purely ingoing waves as r — rp. 
Since for a > M the potential V is regular at any r — 
(including also J'o/Af = 0), we can write 

Vir)^Viro)+0{r-ro). (8) 

By expanding Eq. ([5]) in series around r = ro we find 
that the general solution is a superposition of ingoing 
and outgoing waves: 

r ^ Ae-*'^''^•+Be*"'■•+0(r.-r,(ro))^ fc^ = F(ro) , 

(9) 

where the sign of k is chosen to recover the well-known 
boundary condition for a wave-function in an extreme 
Kerr background (a — s- M and tq — )• M): k ~ uj — rjift, 
where f2 = 1/(2M) is the angular velocity of an extreme 
Kerr black hole. Purely ingoing boundary conditions at 
the stringy horizon imply _B = in Eq. ([9]) or, equiva- 
lently, 

dY 

^-^ikY, r^ro. (10) 

This is the condition we impose in our numerical code. 

For each w, we integrate Eq. ([5]) numerically inward, 
starting at some large radius (typically r^o = 400 Af) 
where we impose the asymptotic behavior ([7]). Our re- 
sults are robust to variations of in a reasonable range. 



We stop the numerical integration at r = rg, where the 
value of the field Y{uj, tq) is extracted. Finally, we repeat 
the integration for different values of uj until the desired 
boundary condition (either F(a;,ro) = or Eq. ([TU)) ) at 
r = ro is satisfied, typically to within an accuracy of 
10-1". 

In our numerical computations we make use of the 
series expansion (j6]), truncated at fourth order. When 
I aw I < 1, the series expansion is a very good approxima- 
tion of the exact eigenvalues. However, in some cases (i.e. 
when I aw I ^1), instead of the series expansion we have 
used exact numerical values of sAim obtained by solving 
Eq. ([3]) with the continued fraction method [l^ . 

We focus on the most relevant gravitational perturba- 
tions, described by the Teukolsky equation with s = 2. 
To compute unstable modes we also make use of the sym- 
metry [ii [Hi 

m ^ -m , w -w* , sAijn s A*-m • (11) 

In practice, this symmetry means that modes with az- 
imuthal number —m can be obtained from those with 
azimuthal number m by changing the sign of the real 
part of the frequency. Therefore we focus on modes with 
Re[w] = w_R > only. 



II. PERFECT MIRROR 

A. Unstable modes with I — m = 2 

Let us start by reviewing and extending the results of 
Ref. 0, which first found that superspinars with a per- 
fectly reflecting surface are unstable due to the ergoregion 
instability. In the top panel of Fig. [T]we show unstable 
frequencies for s = / = m = 2asa function of the spin pa- 
rameter a/M for selected values of ro/A/. We see that the 
instability (signalled by a positive imaginary part w/ for 
the frequency) is always strong, i.e. it always occurs on a 
short timescale r = l/w/ lOA/ - 5 x IQ-^{M/Mq) s, 
at least when a < 2.2M . It is interesting to note that the 
instability is also effective for = M and a = M + e, i.e. 
for an object as compact as an extremal Kerr BH with 
a rotation parameter that only slightly violates the Kerr 
bound. This is also illustrated in Table HI 



From Fig. [T] and Table U it is clear that when ro > M 
the imaginary part does not vanish as a — > M . This is in 
agreement with our expectations, since when ro > rn the 
"BH bomb" instability 0, [Uj occurs even when a < M. 



The dependence of the eigenfrequencies on the mirror 
location is also shown in the bottom panel of Fig. [1] for 
different values of the spin parameter. The imaginary 
part of the frequency is positive (i.e. the object is un- 
stable) for a wide range of parameters. For any value 
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FIG. 1: Top: Real (left) and imaginary part (right) of unstable gravitational modes of a superspinar as a function of the spin 
parameter, a/M, for Z = m = 2 and several fixed values of ro. Bottom: Real (left) and imaginary part (right) of unstable 
gravitational modes of a superspinar as a function of the mirror location, ro/M, for I = m = 2 and different fixed values of the 
spin parameter. Large dots indicate purely imaginary modes. 



TABLE I: Unstable gravitational (s = 2) frequencies with I — 7Ji — 2 for a superspinar with a perfect reflecting surface {TZ = 1) 
and with a "stringy event horizon" (TZ = 0) at r = ro. All modes in this table have been computed using numerical values of 
sAim obtained via the continued fraction method [2^ . 



ro/M 
0.01 
0.1 
0.5 
0.8 
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{ujrM ,ojiM), 7^ = 1 



{ujrM,ujiM), 7^ = 



a = I.IM 
(0.5690,0.1085) 
(0.5548,0.1237) 
(0.4571,0.1941) 
(0.3081,0.2617) 
(0.1364,0.3095) 
(0.0286,0.3248) 



a = l.OlAf 
(0.9744, ,0.0431) 



(0.9673 . 
(0.9256 , 
(0.8598 , 
(0.6910. 
(0.4831 



0.0475) 
0.0631) 
0.0878) 
0.1742) 
0.2655) 



a = l.OOlAf 
(0.9810 , 0.0097) 
(0.9794 , 0.0110) 
(0.9688,0.0155) 
(0.9507 , 0.0202) 
(0.9003 , 0.0640) 
(0.6071 ,0.2207) 



a = I.IM 
(0.5002,0.0173) 
(0.4878,0.0260) 
(0.3959,0.0719) 
(0.2537,0.1053) 
(0.0916,0.1219) 
(-0.0078,0.1233) 



I.OIM 
(0.9498 , 0.0062) 
(0.9435 , 0.0093) 



(0.9016 
(0.8298 
(0.6530 
(0.4377 



0.0237) 
0.0376) 
0.0821) 
0.1230) 



a = l.OOlAf 
(1.0286,0.0033) 
(1.0252,0.0048) 
(1.0052. ,0.0091) 
(0.9793,0.0095) 
(0.8853,0.0313) 
(0.5696,0.1064) 



of a/M in the bottom panel of Fig. [T] the instability is 
strongest when r^/M ^ 1, and is effective also in the 
limit Vq/M <C 1 (although in this regime high-energy 
corrections to the background metric could be relevant). 

Overall, Fig. [T]shovirs that the strongest instability oc- 
curs roughly when a/AI ^ 1.1. For larger values of the 
spin the imaginary part decreases and eventually it van- 
ishes (causing the instability to disappear) for a critical 
value of a/M which depends on vq. At first sight, this re- 
sult seems in contrast with the superradiant nature of the 



instability, as one may naively think that the instability 
should become stronger for large spins. 

In Fig. [5] we show that this expectation is not justi- 
fied by plotting the proper volume of the ergoregion as a 
function of a/M. The proper volume can be computed 
via 

V = 47rJ d9 J dry/grrgeeQw > (12) 
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FIG. 2: Proper volume of the ergoregion as a function of the 
spin a/M. The volume increases monotonically when a < M, 
is infinite at a = M and decreases monotonically when a> M. 
The proper volumes for a ~ 2M and a ~ 0.3M are roughly 
the same. In the inset we show the azimuthal section of the 
ergoregion for selected values of the spin. These spins are 
marked by filled circles and capital Latin letters in the main 
plot; their numerical value is indicated in parentheses in the 
figure. In the limit a/M — > oo the ergoregion becomes so 
oblate that its proper volume shrinks to zero. 



where we have considered a constant time shce, the met- 
ric elements are taken from Eq. ([1]), we have exploited 
the reflection symmetry of the Kerr metric, and we have 
already integrated out the dependence. For a < M 
the ergoregion extends between the outer Kerr horizon 
at th = + V — a? and the "outer ergosphere ra- 
dius" at ri,+ {0) = AI + V — a? cos^ 9. In this case 



we set Ti ~ rn, ff 



' e+, 



and 9i 



in the inte- 



gral above. A straightforward calculation shows that 
in this case the proper volume increases monotonically 
with a/M, eventually diverging^ for a = M. However, 
when a > M , the ergoregion extends between the in- 
ner ergosphere at re~{9) = M — \J A/^ — cos^ 9 and 
the outer ergosphere re+(0). Therefore we set ri = re_, 
Tf = re+, and 9i = arccos(M/a) in the integral ([T^. 
In this case, the proper volume of the ergosphere mono- 
tonically decreases with a/M. In the inset of Fig. [5] we 
plot an azimuthal section of the ergoregion for selected 
values of the spin parameter. The proper volume of the 
ergoregion vanishes as a/M —> oo because the ergoregion 
becomes more and more oblate (in the equatorial direc- 
tion) as the spin increases. As a/M — )■ oo the proper 
volume shrinks to zero and the ergoregion instability for 
modes with I = m = 2 becomes harmless. In Section HVl 
however, we will see that this suppression of the ergore- 
gion instability is less effective for modes with I = m ^ 2, 
which are more concentrated in the equatorial region and 



which make superspinars unstable even for larger values 
of the spin. 



B. Unstable modes with m = 

Supcrradiance due to an ergoregion is not the only 
mechanism driving instabilities in superspinars. Unsta- 
ble modes also exist for m = 0, when the condition for 
supcrradiance uj < mfl = cannot be fulfilled. This is 
shown in the left panel of Fig. [31 where we show different 
gravitationally unstable modes for I = 2 and m = 0, 1, 2. 

We see that the unstable mode with m = 2 exists in 
the range < r^/M < 2, i.e. out to the outer location of 
the ergoregion in the equatorial plane. Unstable modes 
with m ~ 1 and m ~ exist only in a more limited range 
around rg ^ M. Modes with larger values of m(< /) 
drive stronger instabilities, but the instability for modes 
with m = is also important, because it occurs on a 
dynamical timescale r = l/w/ ~ M for a wide range of 
parameters. In the right panel of Fig. [3] we show some 
unstable modes with to = for different values of the 
spin parameter. 

Unstable modes with m = have been recently found 
by Dotti and Gleiser [13, Ell- By imposing regularity con- 
ditions at r = — oo these authors found (an infinite num- 
ber of) unstable purely imaginary modes when a > M. 
The supcrspinar model we are discussing reduces to the 
spacetimc considered in Refs. jTol[T]| when tq — > — oo and 
TZ = I. We carried out a search of these purely imagi- 
nary unstable modes, and our results agree very well with 
those of Ref. [ll| in this limit. For illustration, in the left 
panel of Fig. |3] we show that the frequency of the m = 
purely i mag inary mode for a = 1.4M matches the result 



of Ref. [Ill for ro — > — oo. The figure shows that the 
frequency of these modes settles to its asymptotic value 
when ro < — 3M, so it makes sense to fit Mw/ for these 
m — purely imaginary modes as a function of a, setting 
To = — 3Af. A comparison between the numerical results 
and the polynomial fit 



Muji = 6.375 + 0.177a/M + 0.230ia/Mf 



(13) 



^ This is because when a = M, grr ~ l/{r — M)^ near the horizon 
r = M: this causes the integral II12I I to diverge logarithmically. 



is presented in the right panel of Fig. |4l In the range 
1.15 < a/M < 2 the fit is accurate to within 0.003% and 
suggests that the imaginary part of the frequency (and 
therefore the "strength" of the instability) grows approx- 
imately as a quadratic function of the spin. We stress 
that these results have been obtained using the asymp- 
totic expansion of the angular spheroidal eigenvalues |2l| 

,Ai^ ^{2l-r]+ l)\auj\ - s(s + 1) + O(|aw|0) , (14) 

where 77 = 2 max (|to|, |s|). Because for these modes 
I aw I > 5, this expansion is a good approximation of the 
numerical ei genv alues computed by the continued frac- 
tion method |2li . 
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FIG. 3: Left: Imaginary part of unstable gravitational modes of a superspinar as a function of the mirror location, ro/M, for 
a = I.IM, 1 = 2 and m = 0, 1, 2. Right: Imaginary part of unstable gravitational modes of a superspinar as a function of the 
mirror location, ro/M, for I = 2, m = Q and several values of the spin parameter, o. 
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FIG. 4: Left: Purely imaginary unstable mode as a function of the mirror location, ro/M < 0, for a = 1.4M, s = I = 2 and 
m = 0. In the limit ro — >■ — oo, Mujj ~ 7.07, which perfectly agrees with results in Ref. Right: Purely imaginary unstable 
mode as a function of the spin, a/M, for ro = — 3M, s = I = 2, m = 0. Numerical results (black straight line) are consistent 
with the quadratic fit of Eq. (|13l) (red dashed line). The dot marks the case considered in the left panel. 



III. ABSORBING BOUNDARY CONDITIONS 
(HORIZON-LIKE SURFACE AT r = m) 

From the results discussed in the previous section we 
conclude that a dynamical instability is almost unavoid- 
able in a broad region of the parameter space if the sur- 
face of the superspinar is perfectly reflecting. The in- 
stability is present even in what would naively seem the 
most phenomcnologically viable case, i.e. when tq M 
and a = M + e. One could argue that a perfectly re- 
flecting surface maximizes the efficiency of the ergoregion 
instability because negative-energy modes, which are po- 
tentially dangerous, cannot be absorbed, and that this 
might not happen for different boundary conditions. In 
fact, Kerr BHs are stable because (despite supcrradiant 
scattering) the negative-energy modes can flow down the 



horizon. Therefore we expect ingoing boundary condi- 
tions (7?. = 0) at r = ro to represent the worst possi- 
ble situation for the ergoregion instability to develop. If 
we find an instability even in this case, the superspinars 
described by our simple model are doomed to be unsta- 
ble. This choice also seems more physically motivated 
than the perfectly rcfiecting boundary conditions, be- 
cause ro might be the location of an event horizon formed 
by string-inspired modifications of gravity at high curva- 
tures. We consider both modes with I ^ rn ~ 2 (which 
we expect to be affected by the ergoregion instability) as 
well as modes with m — Q, which we found to be unstable 
in the perfectly reflecting case. 

The punchline of this section is that unstable modes 
exist even when we impose ingoing boundary conditions. 
Qualitatively, the results are the same as those obtained 
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FIG. 5: Top: Real (left) and imaginary part (right) of unstable gravitational modes of a superspinar as a function of the spin 
parameter, a/A/, for I = m = 2 and several fixed values of the horizon location ro/M. Bottom: Real (left) and imaginary part 
(right) of unstable gravitational modes of a superspinar as a function of the horizon location, for I = m = 2 and fixed values of 
the spin parameter. Large dots indicate purely imaginary modes. 



by imposing perfect reflection at the surface of the su- 
perspinar. The instability is slightly weaker than in the 
previous case, but it is again unavoidable in a wide region 
of parameter space. 

Gravitationally unstable modes with I ~ m = 2 and 
TZ ^ are listed in Table |T] and shown in Fig. [S] (to be 
compared with Fig. [1]). Typically the imaginary part of 
the unstable modes when 7?, = is only one order of mag- 
nitude smaller than that obtained imposing TZ ~ 1, which 
causes the instability to disappear at slightly smaller 
spins, i.e. when a/M > 1.75. However, as already men- 
tioned, in the next section we will present evidence that 
higher-/ modes are unstable for larger values of the spin 
and show that our results arc sufRcicnt to rule out su- 
pcrspinars as astrophysically viable alternatives to Kerr 
BHs. 

Also, we stress that perfectly absorbing "stringy hori- 
zons" can only be created by high-energy effects tak- 
ing place beyond the range of validity of general rela- 
tivity. From a phenomenological point of view, the re- 
gion that should be modified by these high-energy cor- 
rections is close to the curvature singularity of the Kerr 



metric {tq <C M). Unstable modes generically exists for 
a = M + e, even in the limit ro/M — 0. Moreover, re- 
sults are smooth in the limit vq/AI — > 0, which means 
that in the region spanned by our calculations curvature 
singularities do not affect our conclusions. 

Finally, Fig. [5] (to be compared with Fig. [3]) shows 
that unstable modes with m = are still present when 
we impose 7?, = at r = tq . 



IV. MODES WITH Z = m > 1 AND PHYSICAL 
ORIGIN OF THE INSTABILITY 

We have seen that when we impose perfectly reflecting 
boundary conditions {TZ = 1) superspinars are plagued 
by several instabilities. Moreover, these instabilities are 
still present when we impose TZ = 0, i.e. when we 
consider a "stringy horizon" at r = ro, which would 
be expected to quench the instabilities. Here we ana- 
lyze in more detail how these different instabilities arise. 
We focus first on the case tq/M > 0. For TZ = 1 
from Figs. [1] and [3] (and analogously for 7?. = from 
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FIG. 6; Imaginary part of unstable gravitational modes of a 
superspinar as the spin parameter, a, for Z = 2, m = and 
several values of the horizon location, tq. 



G 1.0 




FIG. 7: Critical radius rait/M as a function of a/M for s — 
I — 2 and m = 0, 2. We impose both TZ = and TZ — 1 
at r = To. For m = 2 an instability occurs when vq < Va-it- 
For m = an instability occurs in the region delimited by the 
curves on the left. When i = m 3> 1 an instability is expected 
to occur also when a > 2.2M in the TZ = 1 case, and when 
a > 1.75M in the 7^ = case (cf Fig. [8|. Although not 
shown, the region where ro /M < is expected to be unstable 
for any value of a > M. 



Figs. [S] and (5]) we see that, for both m = and m = 
2, the imaginary part vanishes at some critical radius: 
a;/(rcrit) = 0. By using a root-finding routine wc can 
solve for the critical radius as a function of the spin 
parameter a/M . The results are shown in Fig. [71 For 
TO = 2 and TZ = 1 the instability occurs in the region 
below the dot-dashed blue line extending from a ^ M 
up to a/M ~ 2.2. The dashed horizontal line marks the 
location of the outer ergoregion on the equatorial plane 
(r = 2M). As a/M — ^ 1 the critical radius roughly coin- 
cides with the location of the ergoregion, and it decreases 
monotonically for larger rotation. The situation is sim- 
ilar for m = 2 and 7^ = (solid black line extending 



from a ~ M up to a/M ^ 1.75), with the instability 
disappearing earlier (at a/M ~ 1.75) because the ingo- 
ing boundary conditions allow negative energy modes to 
flow down the stringy horizon. This plot confirms our 
qualitative understanding of the instability: as shown in 
Fig. [21 in the limit a/M — >■ oo the proper volume of the 
ergoregion vanishes and superradiance cannot destabilize 
the modes at arbitrarily large spins. 



10' 



10 





i i i 1 i i i 1 i i i 1 
R=1 ^ s ^ J=m=5 


/ — -.4=111=2 \ 


\l=m=4 \ 
\l=m=3 \ \ 




l=m=2 \ \ \ ■ 

\ ^ ^ 


■ R=0 \. \\ 


1 \ t 
I 1 * 
' i ' 




! I . 
1 . i 


l=m=5\ \ ' 


1 ' ' 
. . 1 . . 1 ; . . . I- 



1.0 1.5 2.0 2.5 3.0 3.5 

a/M 



FIG. 8: Imaginary part of the fundamental unstable mode as 
a function of the spin a/M for a superspinar with ro = 0. 
Upper curves refer to 7?. = 1 with Z = m = 2, 3, 4, 5. Lower 
curves refer to 7?. = with Z = m = 2, 3, 4, 5. Dots indicate 
purely imaginary modes. 

It would therefore seem that superspinars with a > 
2.2M might be stable, for any value of tq and even in the 
most restrictive case with 7?, = 1. We argue, however, 
that this is not true if we consider modes with I = m > 
2. In fact the angular distribution of modes with higher 
Z = TO is more concentrated around the equatorial plane. 
Higher-/ modes become more effective at destabilizing 
the superspinar when the ergoregion is oblate, i.e. when 
a ^ M (see again the inset of Fig. [SJ. This intuitive 
understanding is confirmed by Fig. [SI There we plot the 
imaginary part of the fundamental unstable modes for 
/ = TO = 2 and I ~ m > 2 (setting rg = 0), both for 
7^ = 1 and 7^ = 0. 

In the case 7?. = 1, unstable modes with higher I = m 
exist for larger values of the spin. For example the 
I = m ^ 5 mode becomes stable when a/M > 3.6, while 
the / = TO = 2 mode becomes stable when a/M > 2.2, 
as previously discussed. Modes with I = m ^ 1 are 
generally difficult to compute with our code. Our re- 
sults suggest that, for any fixed value oi a/M ^ 1, there 
are always unstable modes as long as / = to is sufficiently 
large. We stress that these results are in contrast with the 
case of Kerr BHs, where the superradiant amplification is 
always stronger for I = m = 2 |22|. Results are qualita- 
tively similar for 7?, = 0, in which case the I = m = i in- 
stability disappears when a/M > 2, while the Z = to = 2 
instability disappears when a/M > 1.75, as previously 
discussed. However in this case the / = to = 4 insta- 
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bility is weaker than the Z = m = 3 instabihty, and the 
I = m = 5 instabihty disappears for smaUer values of 
a/AI than in the I = m = 2 case. This is because the 
stringy event horizon quenches the instabihty of higher-Z 
modes, similarly to the horizon of a Kerr BH. 

More in general, the ergoregion instability that we have 
found here at the linear level can be related to simple 
kinematical properties of the Kerr spacetime. In fact, 
as we show in Appendix [Xj Kerr spacetimes with a > M 
admit stable non-equatorial null circular orbits with nega- 
tive energy. These orbits exist at any radius r < M . The 
very existence of these orbits is enough to prove that the 
spacetime is plagued by the ergoregion instability, pro- 
vided that purely reflecting boundary conditions (TZ = 1) 
are imposed at < M. This is because null orbits are 
the geometric-optics limit of gravitational perturbations, 
as can easily be seen by expanding their propagation 
equation in powers of 1/A (A being the wavelength of the 
perturbation) . Therefore, the existence of stable null cir- 
cular orbits with negative energies implies the existence 
of short-wavelength modes with negative energies. Under 
perfectly reflecting boundary conditions, these modes can 
only leak to infinity by tunneling through the potential 
barrier. However, because particles outside the ergore- 
gion must have positive energies, this leak makes the en- 
ergy of the perturbations inside the ergoregion more and 
more negative. As a result, their amplitude grows with- 
out bound, thus revealing the instability of the spacetime. 
For these reasons we expect modes with Z = m ^ 1 to 
be unstable for arbitrarily large values of a/M, at least 
for TZ = I. This expectation is consistent with the gen- 
eral theorem of Ref . Q , which states that any spacetime 
possessing an ergoregion, but not an event horizon, is vul- 
nerable to the ergoregion instability. As shown in Fig. El 
this expectation is not always justified in the less efficient 
case with TZ = 0. 

The existence of stable non-equatorial null circular or- 
bits with negative energy clarifies why superspinars are 
unstable even under perfectly absorbing boundary con- 
ditions (i.e., in the presence of a stringy horizon), while 
for Kerr BHs with a < M the presence of the horizon 
kills the ergoregion instability. For superspinars, the ef- 
fective potential for gravitational perturbations presents 
a minimum at small radii (corresponding, in the eikonal 
limit, to the location of a negative-energy stable non- 
equatorial null circular orbit), and then rises as r/M ^ 0. 
Therefore, the ergoregion modes need to tunnel through 
a potential barrier to fall into the "stringy" horizon. The 
stability of the supcrspinar depends on a delicate balance 
between the transmission coefficients through the "inner" 
and "outer" potential barriers. 

A possible objection against instability could be the 
following. For sufficiently fast rotation (perhaps even for 
spins as low as a ~ 6M, i.e. at the higher end of the vi- 
able range identified by Ref. [i3|)j if unstable modes exist 
in the eikonal limit [l — m oo) their imaginary part 
will be small, even in the case TZ = 1. The ergoregion 
instability is due to ergoregion modes leaking to infinity 



through the potential barrier, but the tunneling becomes 
less and less effective as the modes behave more and more 
like particles, because the amplitude transmitted to in- 
finity scales as exp(— L/A) (where L ~ M is the width 
of the barrier and A is the mode's wavelength). It is 
therefore conceivable that the imaginary part might be 
so tiny that these modes can be considered stable for all 
practical purposes. 

However, accretion is known to spin superspinars down 
. According to Ref. [2^ , a BH which is initially non- 
rotating gets spun up to the extremal limit a = M, where 
it cannot be spun up any more [26l . [27j , by accreting a 
mass AM = (V6 - l)A/in = 1.4495Mi„ (Mjn being the 
initial BH mass). This corresponds to the accretion of 
a gas mass AMq = 1.8464Min, of which AM falls into 
the BH and AMq — AM is dissipated by the disk's vis- 
cosity into electromagnetic radiation. Similarly, a supcr- 
spinar with a/M = 7 gets spun down to a/M = 1.5M 
(where the ergoregion instability is always effective) by 
accreting a mass AAI = 1.730Min (corresponding to a 
gas mass AMq = 2.295A'/in). The two processes (spin-up 
of a Schwarzschild BH to the extremal Kerr limit, and 
spin-down of a supcrspinar from a/M — 7 to a/M =1.5) 
involve amounts of accreted material of the same order of 
magnitude, hence the corresponding timescales too will 
be comparable. Supermassive BHs are expected to be 
spun up to the extremal Kerr limit by coherent accre- 
tion'^ on a timescale much smaller than the Hubble time 
[3^ , so a supcrspinar should be spun down to the unsta- 
ble region on a timescale much smaller than the Hubble 
time. For this reason, the existence of supermassive su- 
perspinars is unlikely in the real Universe. 

The situation is slightly different for stellar-mass su- 
perspinars. Analytical arguments [37[ and population 
synthesis calculations [s^ show that BHs in binaries es- 
sentially retain the spin they had at birth, so it is unclear 
whether accretion would be efficient enough to destabi- 
lize a supcrspinar. On the other hand, as far as we know, 
no realistic collapse scenario leading to the formation of 
stellar-mass superspinars has been proposed so far. Typ- 
ical equations of state lead to compact stars rotating with 
a/M < 0.7 0, Q. Polytropic differentially rotating stars 
with a « I.IM can in principle exist 0, Q, but they are 
stable. Even if depleted of 99% of their pressure and 
induced to collapse, these stars do not form a BH and 



^ It has been proposed that supermassive BHs may accrete small 
lumps of material with essentially random orientations of the or- 
bital angular momentum. This "chaotic accretion" results (on 
average) in a spin-down of the BH [28l | , so it is very hard to pro- 
duce fast spinning BHs at all (whereas spin estimates as large as 
a+0.009 



0.989:l:0 002 have been reported [29||). Therefore it should 
be even harder to produce superspinars by chaotic accretion. Bi- 
nary BH mergers are also known to always produce spins below 
the Kerr limit [30l435| . so one would be left only with the possi- 
bility of postulating that supermassive superspinars are born in 
the early Universe due to high-energy physics effects beyond the 
realm of classical general relativity. 
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produce cither a supermassive star (which will collapse 
to a BH with a < M when enough angular momentum 
has been shed in gravitational waves) or a stable, rapidly 
rotating star. 

For tq/M > there is a second family of unstable 
modes with m ~ that cannot be superradiant modes. 
Figs. [3] and [H] show that, for fixed values of a/M, these 
modes only exist in a limited range of r^/M . This range 
corresponds to the blue dot-dashed line (7?. = 1 case) and 
to the solid black line (7?. = case) on the left of Fig. [71 
showing that this family of unstable modes only exists 
for a/M < 1.12. 

These unstable modes are related to the existence of 
stable "polar" null circular orbits, i.e. circular non- 
equatorial orbits with vanishing azimuthal component of 
the orbital angular momentum {Lz = 0) [s^] . Eq. (|A22[) 
of Appendix |X] (see also Fig. [9]) shows that for the 
Kerr spacctimc such orbits exist when 1 < a/i\/ < 
(-3-h2A/3) « 1.17996. For Z = 4 and 7^ = 1 the in- 
stability range for modes with ?n = is 1 < a/M < 1.14. 
However the upper limit of this range is a slowly increas- 
ing function of Z, and it is plausible that in the eikonal 
limit it should tend to (-3 + 2^3) 1.18. 

In conclusion, let us discuss the case tq/M < (with 
TZ — 1), summarized in Fig. |4l Now the ring singularity 
at r/M = is naked, and the spacctime also possesses 
closed timelike curves Therefore it is not surpris- 

ing that an infinite number of unstable modes exist also 
at the linear level [ll|. At variance with the crgoregion 
instability, in the present case the imaginary part of the 
frequency (and therefore the "strength" of the instabil- 
ity) grows roughly quadratically with a/M (cf. the right 
panel of Fig. |4]) . The same kind of instability has been 
found in charged, spherically symmetric BHs with naked 
singularities [lOl and therefore it is not related to rota- 
tion, but to causality violation (see also the discussion at 
the end of Ref. [l3|)- As a matter of fact, we could not 
find any mode belonging to this family when tq/M > 0, 
i.e. when the naked singularity is covered. 

In summary: superspinars arc plagued by several in- 
stabilities for both perfectly reflecting boundary condi- 
tions {TZ = 1 at r = j'o) and perfectly absorbing bound- 
ary conditions {TZ ~ at r = rg). The instability of 
modes with I — m is related to superradiant scattering. 
When ro ~ M, unstable modes with m = also exist 
below some critical rotation parameter: this instability 
is related to the existence of stable polar null circular or- 
bits in the spacctimc (cf. Appendix Finally, when 
tq/M < 0, a third family of m = modes exists pT| . 
This third family of unstable modes is probably related 
to the existence of naked singularity and closed timelike 
curves in the spacetime. 



V. CONCLUSIONS 

The results reported in this paper indicate that su- 
perspinars are unstable independently of the boundary 



conditions imposed at the "excision radius" tq and in 
a significant region of the two-dimensional parameter 
space {a/M, tq/M), if not in the whole parameter space. 
The most effective instability at low rotation rates cor- 
responds to the I = m = 2 (superradiant) mode, but 
when a M and rp ^ M unstable modes with m = 
also exist. The I = m = 2 mode eventually becomes 
stable at large rotation rates, but unstable modes with 
I — m ^ 1 are expected to exist for any value of a/M, 
at least for TZ = 1. While the instability timescale of 
higher-/ modes may turn out to be very long, making 
them marginally stable for practical purposes, the low- 
l instability (which affects superspinars with a/M < 2) 
takes place on a dynamical timescale. Accretion is known 
to spin superspinars down [l^ . so our results indicate 
that superspinars are unlikely astrophysical alternatives 
to Kerr BHs. 

One possible objection is that, in order to assume in- 
going boundary conditions at the surface of the super- 
spinar, we must assume that general relativity is modi- 
fied in that region. Such a modification of general rela- 
tivity in the excised, high-curvature region surrounding 
the singularity is implicit in the original superspinar pro- 
posal by Gimon and Horava [ij , who invoke string theory 
in order to violate the Kerr bound a < M. We stress, 
however, that our results hold for a wide class of theories 
of gravity. Many proposed alternative theories of gr avity 
admit the Kerr spacctime as an exact solution [4l|, \^ . 
Among these theories, we focus on the large class consist- 
ing of Brans-Dickc gravity (with or without a potential) , 
and theories that can be reduced to Brans-Dicke the- 
ory with a potential via a conformal transformation (e.g. 
/(i?) gravity, both in the metric and Palatini formal- 
ism 0,143]). All of these theories admit Kerr-(anti) de 
Sitter as an exact solution if the scalar field is constant. 
When perturbed, these solutions satisfy different equa- 
tions in general relativity and in modified gravity theo- 
ries [i^ |45| , due to the presence of an extra scalar degree 
of freedom (the Brans-Dicke scalar) , so one might naively 
expect the stability properties of the Kerr spacetime to 
be different. However, one can redefine the tensor modes 
via a conformal transformation so that the vacuum ten- 
sor and scalar perturbation equations decouple at linear 
order [45|. Basically this happens because the Brans- 
Dicke action reduces to the Einstcin-Hilbert action plus 
a minimally coupled scalar field in the Einstein frame, if 
no matter fields are present Therefore, the tensor 

modes satisfy the same equations in general relativity as 
in Brans-Dicke theory (or in any other theory that can 
be recast in Brans-Dicke form via a conformal transfor- 
mation.) This means that Eqs. ^ and ([3]), which arc the 
starting points of our analysis, retain their validity, and 
therefore that the instability operates in a wider class of 
gravity theories. 

Our results do not imply that there cannot be stable 
ultracompact objects with J/AP > 1. However, they do 
imply that either (i) Einstein's gravity should be modified 
in such a way as to retain Kerr as a solution, while at 
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the same time allowing the tensor modes and the "extra" 
modes to couple at linear order, or (ii) the structure of 
astrophysical superspinning objects, if they exist at all, 
is not described by the simple Kerr-based superspinar 
proposal of Ref. [l| . 
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Appendix A: Geodesies in D-dimensional Kerr 
spacetimes 

The main goal of this appendix is to study the exis- 
tence of stable null circular orbits (SNCOs) with nega- 
tive energies in Kerr spacetimes. We are interested in 
these orbits because the very existence of SNCOs (or 
more generally, the existence of stable null bound orbits) 
with negative energies is enough to show that a space- 
time is subject to the ergoregion instability, provided that 
purely reflecting boundary conditions are imposed at the 
excision surface r = rp. For completeness we consider 
Z?-dimensional Kerr spacetimes with only one nonzero 
angular momentum parameter, and we specialize to the 
"ordinary" D = 4 case at the end. Our main results 
for four-dimensional Kerr spacetimes are summarized in 
Fig. ini The meaning of the different curves on this plot 
is explained below. 

The metric of a D-dimensional Kerr BH with only one 
nonzero angular momentum parameter is given in Boyer- 
Lindquist-typc coordinates by [46j 



ds 
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where E = r'^+a^ cos^ i?, A^, = r'^+a^-Mor^-^ 
denotes the metric of the unit [D — 4)-sphere and M o 
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FIG. 9: The sign of ^ = L^/E as a function oi a/M and the 
radius r^jM of null circular orbits. When a > M, ^ — >■ oo 
at Vc — M. Stable circular null orbits exist when a > M for 
rc < rs, and have negative energy when Vc < M. Regions 
marked with NCOs are those where no circular orbits exist. 
Orbits of constant radius ri are stable "polar" null circular 
orbits. 



and a are related to the physical mass M and angular 
momentum J of the spacetime 



167rM , 



D-2 J 



M ' 



(A2) 



with Ad = {2tt)^'^-^^/'^T[{D-1)/2]. The outer horizon is 
defined as the largest real root of rf^+a'^ — Alor^^ = 0. 



1. Equatorial Null Geodesies 

For null geodesies in the equatorial plane (0 = 7r/2) of 
the spacetime (jAip . the radial geodesic equation reads 



-2-2 



1 



(A3) 

where ^ = L^/E and where the dot denotes derivatives 
with respect to the dimensionlcss affinc parameter. Con- 
ditions for circular orbits are Vcs{rc) = V^'g(rc) = 0. The 
condition Vcs{rc) ~ implies 



„-D-3 



Mr 



(A4) 



for direct and retrograde orbits, respectively. For D = A 
the outer horizon is located at r = r/f = A/ + V AI^ — 
and (of course) the Kerr bound implies a/M < 1. The 
condition V^ff {rc) = then leads to three different solu- 
tions: 
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(A5) 
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where 



1/3 



The three solutions are ah real 

phot 



(A6) 



phot 



Orbits with Tc ~ r^j 
(tc = r^_^°^) correspond to unstable direet (retrograde) 
circular orbits lying outside the horizon, whereas rc_ < 
th and therefore this solution does not correspond to 
physical circular orbits. For a > M there is only one real 
solution 



^phot ^ 



1 + cosh 



- cosh 



(A7) 



corresponding to an unstable null circular orbit. This is 
shown in Fig. [5] It is easy to show that the same quali- 
tative results hold also when D > 5. Therefore no stable 
null equatorial circular orbits exist in Z?-dimensional Kerr 
spacetimes with a single angular momentum parameter. 



2. Non-equatorial null geodesies 

Let us now focus on non-equatorial null geodesies in 
D-dimcnsional Kerr spacetimes with a single spin pa- 
rameter. We shall follow and generalize the approach 
discussed in Chandrasekhar's book W^ . 

The separability of the Hamilton- Jacobi equation in 
Kerr spacetime was proved by Carter, who also discov- 
ered an additional constant of motion Q (the "Carter 
constant") besides the energy, the angular momentum 
and the norm of the four- velocity [40l lisj. The same 
procedure can be easily generalized to D-dimensional 
Kerr spacetimes with a single spin parameter, given by 
Eq. (jAl|) . Our basic equations are 

Rir) (A8) 

2 e2 „^„2 , n;-_r„ I /f „^2l^5-n „2„ 
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(A9) 
(AlO) 
(All) 



(r^ + - o^Ad sin^ 6 + a£,[AD - (r^ + a^)] 



where we use a dot to denote derivatives with respect to 
the dimensionless affine parameter, and where ^ = Lz/ E, 
rj — Q/E are two constants of motion. Notice that the 
angular equation (|A9p does not depend on D, and that 
the equations above reduce to Eqs. (190) and (191) of 
Ref. O when L) = 4. 



planar orbit {9 = 9o =constant) are 6(6*0) = = O'(6'o). 
From Eq. (|A9p we see that these conditions are fulfilled 
on the equatorial plane (6*0 = 7r/2) only if ij = 0. For 
^0 7^ I"/ 2 planar orbits exist if 



77 = — cos'* I 



^ = ±a sin" 



(A12) 



For the "plus" branch of the solutions above, the radial 
equation (|A8[) simply becomes f = ±E. These geodesies 
are unbound and describe shear-free null-congruences 
[47I I . The "minus" branch of the solutions describes non- 
equatorial planar orbits. These solutions only exist when 
D = A and a < M. Moreover they always lie inside the 
event horizon, and therefore they do not correspond to 
physical orbits. For these reasons we do not discuss them 
further. 



3. Non-equatorial, circular orbits 

Since no planar, non-equatorial circular orbits exist in 
D— dimensional Kerr spacetimes, let us focus on non- 
equatorial, circular orbits, i.e. orbits with constant ra- 
dius but which are not plan ar (i.e. 9 is not constant). 
These orbits are periodic |49l - [51| and they are often called 
"spherical orbits" in the literature, but here we adopt the 
term "circular orbits" as in Refs. (49l - [5l| . 

The conditions for null circular orbits, R{rc) = and 
R'{r^) = 0, read 

rt + {a'-e-v)rl (A13) 
+ Md[v + (^ - a)']r,'"^ - a^^ = , 

Arl + 2(a^ - - (A14) 
-f (5-i5)A/z3[77+(e-a)>r^ = 0, 

which can be solved for ^ and 77 as functions of re- There 
are two sets of solutions: 



V = 2 ' 



(A15) 



and 



a^{D - 5)MDrt + {D - l)Mcr5 - 2r;P(a2 + r^) 



a{D - b)MDrl + 2arf 
77 = [a{D-h)MDrl + 2ar^Y'^ 



X 



(A16) 



a. Proof of no planar bounded orbits in D dimensional 
Kerr spacetimes 

A relevant question is whether non-equatorial planar 
orbits exist in these spacetimes. The conditions for a 



The first set of solutions implies 9 = constant and indeed 
reduces to the "minus" branch of solutions (jA12[) . which 
do not correspond to physical orbits. 

The second set of solutions, Eqs. (|A16|) . can de- 
scribe bound orbits. The condition of stability is simply 
R"{rc) < 0. By differentiating Eq. (jA8[) twice and using 
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Eqs. (jA16|) wc obtain the following expression for R"{rc): 

R"{r,) ^ [{D-mDrl+2r?]-^ (A17) 
X {-8p-5)(Z?-l)M|,r8 + 32r2(^+i) 

+ 16(1? - 5)A/Drf +^ [a\D - 3) + {D - l)rl] } . 

The stability of null circular orbits depends on the sign 
of the expression above. It is possible to show that stable 
circular orbits exist for _D = 4 and a > M, but not for 
D > 5. Therefore in the following we will specialize to 
D = A spacctime dimensions. 

a. D — 4 Kerr spacetime. 



When D = 4, Eqs. (^161) read 



a{rc - M) 
rl [Aa'^M - r^rc ~ iMf] 
a2(rc -M)2 



(A18) 
(A19) 



These equations correspond to Eqs. (224) and (225) of 
Ref . |47[ , and they can be used to define the shadow cast 
by Kerr BHs or superspinars When = 0, from 

Eq. (|H9|) we have 4a^M - rc(rc - SM)^ = 0, which 
defines the equatorial orbits (|A5P - (|A7|) . In general, how- 
ever, the constant of motion can be positive or negative. 
When 77 < 0, Eqs. (|Il8)) - (|JT9| . together with Eq. (|X9)) 
for the ^-motion, implies that orbits of constant radius 
are not allowed [43]. When 77 > circular orbits are al- 
lowed, and according to Eq. (jAlQp they must satisfy the 



condition Aa^M 
condition reads 



rc{rc -Miy > Q. For a < M this 



phot ^ 



< r 



phot 



where r5^'^°* and r^°*^ refer to retrograde and direct un- 
stable photon orbits in the equatorial plane [Eq. (|A5|) ]. 
More importantly for the analysis of superspinars, when 
a > M the condition 77 > reads 



Tc < r 



phot 



(A20) 



where r^°^ is given by Eq. (jA7p . Notice that the con- 
dition above includes also the singular case Tc = M (in 
fact ^, r/ — > 00 when a > M and j'c M). 

When D = 4, from Eq. (jAlTl we see that stable cir- 
cular orbits exist whenever the orbital radius Vc satisfies 
the relation Tc < rs, with [l^ 



rs_ 
M 



.1/3 



M 



for a> M , 
1/3 (A21) 
for a < M . 



Null circular orbits with radii smaller than this critical 
radius are stable. When a > M, rs < r^°^, i.e. stable 
circular orbits are allowed. When a < M the critical 
radius rs is covered by the horizon, and it becomes "vis- 
ible" to external observers only when a > M. Therefore 
stable null circular orbits may exist only for a > M, while 
orbits with r < rs around BHs with a < M do not have 
a physical meaning because they lie inside the horizon. 

By substituting Eq. ((MT|) into Eqs. (jUSl) and (|JT9)) 
we can compute the corresponding critical parameters 
Virs) = Vs and £,{rs) = ^s- 



Vs 


3M2 






es 


a 


M 






For a given value of a/M, when rj = rjs and ^ = £,s we 
have a marginally stable orbit. If instead ^ we have 
a stable orbit, while rj <r]s gives a stable orbit if a > 3M 
and T] > rjs gives a stable orbit if a < 3M. However, 
these are only sufficient conditions, because other stable 
orbits may exist, far from the critical values 775 and ^g. 
In fact, depending on the value of the spin we can have 
different situations: (i) for a < 3M, if r/ < 775(a) there 
is only one stable circular orbit (with Tc < M), while for 
77 > 77s (a) we have two stable orbits: one with rc < M 
and one with rc > M; (n) for a > 3M, when 77 < 27M^ 
we have only one stable circular orbit (with Tc < M); 
when 27 < 77 < 7^5(0) we have two stable orbits with 
Tc > M and one with rc < M; when 77 > 775(0) we 
have one stable circular orbit with r^ < M and one with 
Tc > M. This can be understood by plotting 77 as a 
function of r, with < r < rs{a), for various values of a. 

Also, let us consider the sign of the impact parameter 
^ — Lz/E. A study of Eq. (jAlSP shows that there is a 
critical spin 



such that: 



a/3 ( -3 + 2V3]m w 1.17996M, 



(A22) 



• if a > fli, then ^ > for rc < M and ^ < for 
rc > M. Notice that ^ diverges if rc = M. 

• if A/ < a < gl, then ^ > for rc < M and for 

ri(a) < rc < r2(a) (with ri,r2 > M), whereas 
^ < elsewhere. Notice that ^ diverges if rc = M. 

• if a < M then ^ > for r+ < rc < r2(a), where r_|_ 
is the outer Kerr horizon and r2 < 3M. ^ < for 

7-c > 7-2(0). 

The situation for a four-dimensional Kerr spacetime is 
summarized in Fig. |9l Orbits of radius ri and r2 carry 
vanishing angular momentum (L^ ~ 0) and therefore are 
called "polar" null orbits. Orbits of constant radius r2 
arc unstable polar null orbits, while orbits of constant 
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radius ri are stable polar null orbits, and they exist for 
M < a < gl. 

A relevant question to ask is whether the null circular 
orbits that we have identified have positive or negative 
energy. The sign of the energy is determined by impos- 
ing that the geodesies be future oriented, i.e. that the 
derivative t of the coordinate time with respect to the 
affine parameter [given by Eq. (jAll|) ] be positive. (This 
is because the hypersurfaces t = const are spacelike for 
any r > if a > M, and for any r > H a < M .) 
By substituting Eq. (|A18|) into Eq. (jAlip . we find that 
for the non-equatorial null circular orbits that we have 
identified we have 



t 



E 



' ' a^cos^i 



M 



(A23) 



Because these orbits cross the equatorial plane (as they 
have ?7 > 0), we can evaluate Eq. (|A23|) for 6 = 7r/2. The 
energy E is a. constant of motion, so it cannot change 
sign along the trajectory. Then it is clear that all orbits 
have > 0, with the exception of orbits with r^/A/ < 1, 
which, as we have seen, exist only for a/M > 1. Indeed, 
it possible to show explicitly that orbits with r^M < 1 
in Kerr spacetimes with a/M > 1 have negative energy 
at all times. Using Eq. (|A9p . one obtains that such orbits 



have polar angle 9 oscillating between tt/2 + 9c and tt/2- 
9r, with 



cos 9r = 



2rc V^-^A {a^M + rl{2rc. ~ MI)) - p 



a^{rc-M)^ 



(A24) 



where p 



SM'^r'^ + 2a'^Mrc. One can show that 



cos^ 6'c < 1 for a/M > 1 and rc/M < 1. Using this 
expression in Eq. (jA23[) it is then possible to show that 
the energy must be negative all along trajectories with 
a/M > 1 and rc/M < 1. The region where stable 
negative-energy orbits exist is shown in Fig. (S) 

Finally, let us suppose we have a compact object ro- 
tating with a > M. According to the cosmic censorship 
conjecture, the singularity at r/M = must be excised 
by some horizon-like one-way membrane or by a reflecting 
surface. It is then natural to ask the question of what the 
compactness of the object can be if one wants to excise 
all SNCOs with negative energies. Because such orbits 
exist for any Tc < A/, if a > M, the maximum allowed 
compactness turns out to be /imax = M/r = 1. Because 
orbits with Vc ^ M lie far away from the singularity at 
r/M = 0, this maximum compactness is not expected to 
be altered by high-energy corrections. 
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